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Abstract 

Mesoscopic density functional theory for inhomogeneous mixtures of sperical particles is devel- 
oped in terms of mesoscopic volume fractions by a systematic coarse-graining procedure starting 
form microscopic theory. Approximate expressions for the correlation functions and for the grand 
potential are obtained for weak ordering on mesoscopic length scales. Stability analysis of the 
disordered phase is performed in mean-field approximation (MF) and beyond. MF shows existence 
of either a spinodal or a A-surface on the volume-fractions - temperature phase diagram. Separa- 
tion into homogeneous phases or formation of inhomogeneous distribution of particles occurs on 
the low-temperature side of the former or the latter surface respectively, depending on both the 
interaction potentials and the size ratios between particles of different species. Beyond MF the 
spinodal surface is shifted, and the instability at the A-surface is suppressed by fluctuations. We 
interpret the A-surface as a borderline between homogeneous and inhomogeneous (containing clus- 
ters or other aggregates) structure of the disordered phase. For two-component systems explicit 
expressions for the MF spinodal and A-surfaces are derived. Examples of interaction potentials of 
simple form are analyzed in some detail, in order to identify conditions leading to inhomogeneous 
structures. 
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I. INTRODUCTION 



One of the main problems that arise in theoretical description of complex fluids is the 
role of density fluctuations on the mesoscopic length scale. Such fluctuations are not im- 
portant in the case of simple fluids, and for this reason simple liquids can be accurately 
described by the liquid theories [l| which focus on the microscopic length scale, whereas the 
long-range fluctuations are treated via mean-field (MF) approximation. An exception is the 
critical region where the long-range fluctuations dominate. Universal features of the critical 
phenomena are described by the phenomenological Landau- Ginzburg- Wilson (LGW) theory 
2|, |3|, because the field-theoretic methods allow for more accurate treatment of the domi- 
nant long-wavelength density fluctuations. In the LGW theory the microscopic structure is 
entirely neglected, however - the pair correlation function in the LGW theory decays mono- 
tonically. The two approaches (i) accurate description of the microscopic length scale and 
rough approximation for the long-wavelength fluctuations, and (ii) accurate description of 
the long-wavelength fluctuations with neglected microscopic structure are complementary. 
Both, used separately, give satisfactory description of simple fluids. Exact theories that are 
capable of description of nonuniversal features of phase transitions were also developed |4t9|, 
but so far these theories were applied to homogeneous phases. 

When there are competing tendencies in the interaction potentials, then self-assembly into 
different aggregates, living polymers, clusters, micelles or another objects ('supermolecules' 
having characteristic size) may occur. A notable example of such interactions is the effec- 
tive short-range attraction long-range repulsion (SALR) potential 10l-ll6|. In addition to 
the liquid order on the microscopic length scale (described by the pair distribution function) 
ordering on the mesoscopic length scale may be present in such systems. This additional 
ordering is associated with packing of the 'supermolecules' in the lyotropic liquid crystalline 
phases. In the liquid theories, designed for description of the microscopic structure, the 
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such additional ordering is manifested by a lack of solutions of the associated 
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15| . On the other hand, in the Landau theory modified by Brazovskii 



the dominant fluctuations of the order parameter (OP) are of finite wavelength 17|. The 
functional of the form postulated by Brazovskii was used for a description of various am- 



phiphilic systems 



18|-|23|, since the dominant finite- wavelength fluctuations of the abstract 



order parameter (OP) can represent in particular the density fluctuations on the mesoscopic 
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length scale. Indeed, when the values of the phenomenological parameters in the functional 
are properly adjusted, the LB theory predicts stability (or metastability) of lyotropic liq- 



uid crystalline phases observed in amphiphilic systems |l8l . |20H23|. One should note that 
in the LB theory fluctuations of the OP lead to a change of the continuous transition be- 
tween disordered and lamellar phases obtained in MF to weakly first order transition [l^ 
which occurs at lower temperature. The long-wavelength fluctuations in the LGW theory 
just modify the critical exponents, whereas the order of the transition remains the same. 
The Landau- Br azovskii (LB) functional is quite general (the OP can have different physical 
meaning), and is expressed in terms of phenomenological parameters (coupling constants) 
whose precise relation to measurable quantities is not known - it should be derived from 
more fundamental microscopic theory. Because the LB theory correctly describes the qual- 
itative properties of systems self-assembling on the mesoscopic length scale, it is desirable 
to find a relation between the Brazovskii theory and the exact statistical mechanics. The 
main questions are: (i) what is the range of validity of the LGW and LB theories (ii) for 
what kind of interaction potentials the LB rather than the LGW theory is valid (iii) how 
the phenomenological parameters are expressed in terms of thermodynamic parameters and 
(effective) interaction potential. The approximate theory derived from the statistical me- 
chanics should allow for determination of phase diagrams and structure when the interaction 
potentials are known. 

A natural theory to start with is the density functional theory (DFT) 2J] which allows 
for description of inhomogeneous systems and in principle is exact. However, the exact form 
of the grand potential functional is not known. In the widely used versions of the DFT the 
contribution to the grand potential associated with interactions is of the MF type. This 
approximation works well for a description of the microscopic structure and away from the 
critical point. However, more accurate approximation for the grand potential functional is 
desirable when the mesoscopic scale fluctuations dominate and may affect the order and 
location of the phase transition to a liquid crystalline phase. On the other hand, we have to 
make simplifying assumptions to make the theory tractable. 

Since we need a description of the ordering on the mesoscopic length scale, we may in- 
troduce mesoscopic density that describes the distribution of particles less accurately than 
the microscopic density, but more accurately than the average density. Such an approach 



was proposed in Ref. 



25l | for a one component system of spherical particles. The meso- 
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scopic density is defined as the microscopic density averaged over regions larger than the 
molecules and smaller than the characteristic length of ordering (for example, the size of the 
clusters). Precise definition is given for multicomponent systems in the next section. Proba- 
bility of spontaneous appearance of particular mesoscopic density field was derived from the 
statistical mechanics by integrating the probability distribution over all microscopic states 
under the constraint of fixed mesoscopic density field under consideration. This method is 
analogous to integrating the probability distribution over all microscopic states under the 
constraint of fixed average density in macroscopic parts of the system. The only difference 
is that the constraint imposed on the microscopic density has a form of the field which 
varies on the mesoscopic length scale. The grand potential functional of the mesoscopic 
density field derived in Ref . I2J] consists of two terms. The first term contains contributions 
from fluctuations on the microscopic length scale under the constraint of fixed mesoscopic 
density p(r). This term resembles standard DPT. The second contribution is associated 
with mesoscopic fluctuations 0(r) that can occur in the system when the constraint p(r) is 
removed. 

In the MP approximation the contribution to the grand potential associated with 
mesoscopic-length scale fluctuations is neglected. In this version of MP the average den- 
sity is approximated by the most probable mesoscopic density 25|]. However, in parts of the 
phase diagram that are close to microphase separation the fluctuation contribution can be 
comparable to the first term, and the average density can be significantly different from the 
most probable mesoscopic density. In such cases the fluctuations cannot be neglected. 

It is worthwhile to note that the relation between the systems with and without the con- 
straint on the mesoscopic density distribution resembles the relation between the canonical 
and the grand canonical ensembles. There is some loose analogy between the system in the 
presence of the mesoscopic constraint imposed on the microscopic density distribution, and 
a macroscopic system with fixed number of particles N = Nq, described by the canonical en- 
semble. When the constraint of compatibility between the microscopic density distribution 
and the mesoscopic field p(r) is removed and (0(r)) = 0, then the system is analogous to 
the open system with fluctuating such that (A^) = A^^o (grand canonical ensemble). The 
canonical and grand canonical ensembles with Nq = (N) are equivalent only far from phase 
transitions, when the fluctuations are small, {{N — Nq)'^) oc xtNq. Close to phase transitions 
the compressibility xt is large (diverges at the transition), and fluctuations cannot be ne- 
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glected. At the phase coexistence the most probable density distribution in the open system 
corresponds to either the gas pg or the hquid pi density in an absence of any constraints 
or external fields. However, in the grand canonical ensemble the ensemble average in an 
absence of any constraints or external fields yields a constant density {pg + pi)/2, although 
homogeneous microscopic states with such density occur with negligible probability. The 
difference between this case and microphase separation concerns the extent of the regions 
having different density (or composition) - mesoscopic rather than macroscopic parts of the 
system - and in turn the time scale associated with displacements of these regions, i.e. with 
the mesoscopic rather than macroscopic fluctuations. While it is justified to neglect macro- 
scopic fiuctuations in studies of coexisting homogeneous phases, in the case of microphase 
separation the mesoscopic fluctuations influence the experimentally observed properties of 
the system. 

The fluctuation contribution to the grand potential reduces to the form similar to the 
LB theory in the case of weak ordering on the length scale signiflcantly larger than the 
molecular size 25[. Such kind of ordering occurs in soft-matter systems, and the results of 
Ref. {25!] conflrm validity of the LB theory for soft matter. The fluctuation contribution can 



25j the DFT and 



be treated by fleld-theoretic methods, and in the theory developed in Ref. 
fleld-theoretic methods are both used. 

The theory developed in Ref. 25| is restricted to a one-component system, whereas the 



soft-matter systems are usually multicomponent. The size of solvent molecules can be several 
orders of magnitude smaller than the size of proteins, nanoparticles or colloids, and the 
solvent molecules can be taken into account only via solvent-mediated effective interactions 
between solute particles 0, l^^]. However, the effectively one-component system might 
lead to incorrect predictions when the size ratio is not very large, and when mesoscopic 
fluctuations of the solvent are important. The purpose of this work is an extension of the 
mesoscopic DFT to the case of multicomponent systems of particles of arbitrary sizes. 

In sec. 2 general framework of the theory for multicomponent systems is introduced by 
a systematic coarse-graining procedure. Mesoscopic volume fractions are deflned, and ex- 
pressions for the grand potential and correlation functions are derived in the same section. 
In sec. 3 approximate theory for weak ordering is developed, and the role of mesoscopic 
fluctuations for stability of the disordered phase is discussed. Two-component systems of 
particles of different sizes are studied in more detail in sec. 4. The MP theory is illustrated by 
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three simple examples. Equations derived in sec. 4 are of MF type, and provide information 
whether inhomogeneities on the mesoscopic length scale (clusters or soft crystals) are formed, 
or the system can phase separate into homogeneous phases for given interaction potentials 
and size ratios. In future studies the theory can be applied to different inhomogeneous 
systems along the lines described in sec. 3. 

II. COARSE GRAINING 

A. microscopic density and microscopic volume fraction 

We consider an n-component mixture of nearly spherical particles, with the components 
labeled by Greek letters. The diameter of the hard core of the particle of the specie a is 
denoted by cTq,. A microscopic state is defined by the positions of the centers of mass of 
Ni, ...,Nn particles, 

-M = {{rr}i=i,...,iv.,« = l,...,n} (1) 

where r" denotes the position of the i-th particle of the a-th specie, and Na denotes the 
number of particles of the a-th specie in the considered microstate. Microscopic density of 
the a-th specie is given by the standard definition, 

p^{r,M):=Y.S{r-vt). (2) 

1=1 

For particles of different sizes, for example for a mixture of nanoparticles and small 
molecules, it is convenient to introduce microscopic density that takes into account distri- 
bution of matter inside the molecules, and instead of ([2]) we introduce 

U^,M):=j:U\r-rt\) (3) 

i=l 

where spherically-symmetric structure of molecules is assumed, J,. fa{r) = Va with Va = 
Ticr^/G denoting the volume of the particle of the specie a, and /a(r) describes distribution 
of matter inside such particle at the distance r from its center. For brevity we shall use the 
notation = J dr, indicating the integration region 5* by J^/^zs when necessary. For constant 
density inside the particle Eq.([3]) reduces to the microscopic volume fraction defined by 

Ur,M):=Y.O{^-\r~rt\) (4) 
i=i ^ 
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where 9{r) is the Heaviside unit step function. Integration of Ca(r,A^) over the system 
volume gives the volume occupied by the particles. The interaction energy for a pair of 
particles i,j of the a, j3 species with the centers at and respectively is 

t^./.(|rr-rjl)= / / /.(|r-rf|)K.,(r,rO/,(|r'-rf|). (5) 

Jr Jr' 

Summation convention for repeated Greek indexes is assumed above and in the whole article. 
V^^(r, r')(ir(ir' is the interaction potential between the infinitesimal volumes dr and dr' 
around the points r and r' inside the particles a and /3. The energy of the system in the 
microstate defined by ([3]) or (j2]) can be written as 

E[M] =^-JJj^{r,M)V^p{r,r%{r\M) (6) 

in n Na 

= ^EEEEf^«^(|r?-rf|) 

^ a=l 13=1 1=1 j=l 



B. Mesoscopic density and mesoscopic volume fraction 

Let us choose the mesoscopic length scale R and consider spheres Sr^t) of radius R and 
centers at r that cover the whole volume V of the system. In order to describe ordering on 
the length scale A, we should choose R < \. We define the mesoscopic density of the specie 



a by an extension of the definition introduced for a one-component system in Ref. 25 1 



P.(r):=T^ / ^^Pa{r\M), (7) 

where Vs = AnR^/S is the volume of the sphere SrIt). Similarly, the mesoscopic volume 
fraction of the specie a at r is defined by 

Ur):=^f Ur',M). (8) 

For an illustration, the mesoscopic density and the mesoscopic volume fraction are shown 
in Fig. 1 for a one-component system, when a single particle is located at r = 0, for three 
different mesoscopic length scales R. In this case the center of the particle is inside (outside) 
the sphere Sr^t) for r < 2R/cra (r > 2R/aa), therefore for r = 2R/aa the number density 
([7]) has a discontinuity (Fig.2). For increasing length scale of coarse-graining, the difference 
between ( and pv decreases. 
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FIG. 1: One-component system and the microstate in which a single hard sphere of a radius (t/2 
is located at r = is considered for different length scales of coarse-graining. The mesoscopic 
volume- fraction defined in Eq.Q is shown by the dashed lines, and pv is shown by the solid 
lines. V = vrcr^/G and p is defined in Eq.Q (based on the standard definition of the microscopic 
density ([2])). Top panel: 2R/a = 1. Central panel: 2R/a = 2. Bottom panel: 2R/a = 5. In this 
particularly simple case the fields ([8]) and ([7]) are functions of the distance r from the center of the 
hard sphere. In each case /^.^(r) = v. The distance r is in a/2 units, C and pv are dimensionless. 
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FIG. 2: Cartoon showing distances corresponding to the discontinuities of the mesoscopic number 
density for different scales of coarse graining. From top to bottom 2R/(t = 1, 2R/a = 2 and 
2R/a = 5, as in Fig.l. Shaded circle represents the particle located at r = 0, open circles represent 
the spheres Sii{r) with the centers at r over which the density or volume fraction is averaged. The 
mesoscopic density at the point r is the number of centers of particles inside the sphere SR{r) 
divided by The mesoscopic volume fraction is the fraction of the volume AttR^ /"i that is 

occupied by the particles. 

For a chosen length scale R the mesostate can be defined by {C} = {Ci(r), Cnij)} or by 
{p} = {Pi(i")! •■•5 Pn(i")}- Note that Eq.(I7]) or Eq.(IHl) describes a constraint imposed on the 
microscopic states, and {C} is equivalent to the constraint ([8]) imposed on all the components. 
For a chosen length scale R all microscopic states can be separated into disjoint subsets, 
such that the microstates belonging to a particular subset are compatible with the same 
constraint (Eq.([8])). Microstates belonging to different subsets are compatible with different 
constraints, i.e. with a different form of {C} or {p}. 

Probability density of a spontaneous occurrence of the mesostate {C} is given by 



(9) 



where 



and 



j DC,... J DQe- 

g-/3(£;[.M]-/^McG(r))_ 



(10) 



(11) 



lMe{{C},R} 

E[J^] is the microscopic Hamiltonian, and JMe{{c},R} ^ symbolic notation for the integra- 
tion over all microstates compatible with {(} according to Eq.Q. fla = fJ^a/va and T are the 



chemical potential of the specie a (in appropriate units) and temperature respectively, and 
/3 = l/ksT, with ks denoting the Boltzmann constant. f2co[{C}] is the grand potential in 
the presence of the constraints {(} (Eq.([H])) imposed on the system. The functional integral 
/ DC^i... j D(n in Eq. fnU]) is over all mesostates {C}, which is indicated by the prime. In 
analogous way we can consider mesoscopic theory based on the mesoscopic density. 

We obtain a mesoscopic theory with the same structure as the standard statistical me- 
chanics. The integration over all microstates is replaced in Eq. ffTOj) by the integration over 
all mesostates. The Hamiltonian is replaced in Eq.(|9]) and ( ITO|) by the grand potential in the 
presence of the constraint of compatibility with the given mesostate that is imposed on the 
microstates. The above formulas are exact. So far we just rearranged the summation over 
microstates. The reason for doing so is the possibility of performing the summation over 
the mesostates and over the microstates compatible with a particular mesostate by different 
methods. 

Grand potential in the presence of the mesoscopic constraint can be written in the form 

n,o = U-TS-p^N^, (12) 

where U, S and Na are the internal energy, entropy and the number of molecules of the specie 
a respectively in the system with the constraint (|H]) imposed on the microscopic densities. 
U is given by the expression 

Um] = l[ I V:°^{r, - r,)UrMr,) = I [ [ U^^ir, - v,)p^{v,)pp{v,), (13) 

where 

V:"^{r, - r,) = V^p{ru)g%"{r, - r^) (14) 

ri2 = ki - r2|, and 



C.(rOC.(r.) ^''^ 
is the microscopic pair correlation function for the volume fraction, in the presence of the 
constraint ([H]) imposed on the microscopic states. ?7^^(ri — r2) is given by an expression 
analogous to Eq.(IT4l). with g"^^ replaced by (7^^, the standard pair correlation function in 
the presence of the constraint ([7]). Note that the above functions differ from each other. 
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In particular, for g^^^ smooth increase from zero is expected for ri2 increasing from zero, 
whereas the correlation function for the microscopic density (Eq. ([2])) vanishes for ri2 < (Jap- 
The advantage of Caij) is its continuity (see Fig.l). In addition, Ca(r) < C^cp for all r, 
where Ccp is the close-packing volume fraction, and the gradient of C, is small, |VC| < 1/R- 
The disadvantage of Ca(r) is the expression for the energy (fT3|) in terms of the pair correlation 
function for the volume fraction, Eq. ffTSj) . which was not studied. The mesoscopic density 
([7]) has discontinuities (see Fig.l), and for significantly different sizes of particles the number 
densities for different components may differ by several orders of magnitude. On the other 
hand, the expression for the energy (1131) has a standard form in terms of the well known 
correlation function. When the ordering occurs on the length scale significantly larger than 
the size of the particles, we can make the approximation (see Fig.l) 

(a{r) ^ pa{r)Va. (16) 

Inserting the above expression for pa into Eq. f|T3l) . yields 

V-,ir,r')^^^^f-^. (17) 



It is important to remember that the approximation f[T71) is only valid when the ordering 
occurs on the length scale significantly larger than the size of particles. 

We further assume that the entropy S satisfies the relation —TS = F^, where is 
the free-energy of the hard-sphere reference system with the constraint (IHl) imposed on the 
microscopic volume fractions. 

C. Grand-potential functional and mesoscopic correlation functions 

Let us introduce external fields {J} = { Ji(r), J„(r)} and the grand-thermodynamic 
potential functional 

Q[{PJ}] := -keT log [ J ^Ci-/ 0(^6'^^^^''^^'^^^- i'^-^''^'^-^'^^]. (18) 

—(3Q[f3J] is the generating functional for the (connected) correlation functions for the meso- 
scopic volume fractions, 

lUAriJ-WlM; 5(/3J.,(rO)...5(/3J«„(r„))- ^ ^ 
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We introduce the notation 



G'^.lri - r^) = (C(ri)C/3(r2))^"" = (Ca(ri)C,(r2)) - {Uri)){Ur,)). (20) 

The relation between the mesoscopic and the microscopic correlation functions resulting 
from the definition of the mesoscopic volume fraction ([8]) is given by 

(C«x(ri)-...-Ca„(r„)) (21) 
_ /' DC,... r DCne-^["-'^^^>'-X^°^'-^^"('-\.,(rO ■ ... ■ C.„(r„) 

= , -T^r , (4.(r')-...-4„(rW)), 

Vs -'r'G5_R(ri) Vs Jr(")GSfl(r„) 

with analogous relation for the correlation functions for the microscopic and the mesoscopic 
densities. Note the difference between the micro- and the mesoscopic correlation functions 
resulting from the integration of the former over mesoscopic volumes. In particular, for 
2R > aai3 the mesoscopic two-point correlation function for densities, analogous to Eq.( !2TI) . 
does not vanish for = Y2. This is because for r' G Sr{yi) and r" G Sr{{yi), such that 
|r' — r"| > (Jq,^, the corresponding microscopic correlation function on the RHS of an equation 
analogous to Eq. fl2T]) does not vanish and contributes to the integral. We should stress that 
in Ref. 25| the theory is based on the mesoscopic density, but for significantly different sizes 
of particles the volume fraction is more convenient, as discussed in the previous subsection. 
Let us introduce the Legendre transform 

m{Q] ■■= miPJ}] + / /3^a(r)a(r) (22) 

where 

is the average field (volume fraction) for given {J}. The equation of state takes the form 

In general C, may differ from any mesostate defined in Eq.([H]). We extend the functional 
VLco beyond the set of the mesostates. Let the extension be defined in Eq. f[T^ on the Hilbert 
space of fields that fulfill the restrictions following from the properties of the mesostates, and 
let us keep the notation fico for this extension. The key restriction on the mesoscopic volume 
fraction is the magnitude and the gradient. In Fourier representation we shall consider the 
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functions that vanish ioi k > n / R, where k is the wave number. As discussed in Ref. 25|, the 
fields with magnitudes exceeding the close packing (such fields belong to the Hilbert space, 
but do not represent any mesostate) are irrelevant, since the corresponding Boltzmann factor 
is very small. We introduce the functional of the form 

/3F[{C}] = /3fi,4{C}] - log [ / / D0„e-^I^/--X^'^«*"('-)11, (25) 



(26) 



where 0a (r) is the local fiuctuation of the volume fraction of the component a, and 

Hfiuciic}, m = ^coiic + 0}] - ficoKc}]- 



We introduced the notation {(f)} = {0i(r), 0„(r)}. From Eqs. fl2^ and ffTSl) it follows that 
the functional fl25l) equals the grand potential, when {(} = {(}, with {(} determined from 
Eq.([24D. By definition (0^) = when {(} = {(}. 

Note that from Eq. fl25l) it follows that the inverse correlation functions (related to the 
direct correlation functions) defined by 

^"^'•-"^ "^"5C«.(ri)...5C..(r„) ^'^^ 

consist of two terms: the first one is the contribution from the fiuctuations on the microscopic 
length scale (< R) with frozen fluctuations on the mesoscopic length scale. This term is 

5"/3^^co[{C}] 



(28) 



(5Cai(ri)...5Ca„(r„)' 

The second term is the contribution from the fluctuations on the mesoscopic length scale 
{> R). From Eqs. (fT8l)-( l27|) we obtain equations relating the inverse correlation functions 
with the many-body correlation functions. In the lowest nontrivial order beyond the mean- 
field approximation and for Jq, = we obtain (see |25| ) 



I I G™„,(ri,r2)C.«.(ri,r2,r) = 0, 

JY^ Jr-2 



(29) 



and 



2C™^(ri,r2)=Q(ri,r2) + 



Eq. ( 129|) is the minimum condition for the grand potential. In the MF approximation 
the second term in Eq. (l29|) is neglected. Since there may exist several local minima, the 
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solution corresponds to a stable or to a metastable phase when the grand potential assumes 
the global or the local minimum respectively. The solution of Eq. fl29|) corresponding to the 
global minimum gives the average density for given n and T in the lowest nontrivial order 
beyond MF. 

In order to obtain the two-point inverse correlation function from Eq.f l30p . Hfiy_c[{C}, {0}] 
in Eq. (!26|) is expanded in (pa, and the expansion is truncated. Since the volume fractions are 
less than unity, the corresponding fluctuations are small and such an expansion is justified. 
In this way an equation relating the two-point inverse correlation function with many-body 
correlation functions is obtained. Approximate equation that can be solved in practice will 
be derived in the next section. From Eqs.f l24p .f l23p and ( I19p we obtain the analog of the 
Ornstein-Zernike equation 

/ CZM^r,)GZ,ir2,r,) = 6{r, - r,)6^^\ (31) 

D. Periodic structures 

Let us consider periodic density profiles 

a(r)=CO + $„(r) (32) 

where 

$„(r + P) = $,(r) (33) 

and P = ^iPi where pj are the vectors connecting the centers of the nearest-neighbor 
unit cells and Ui are integer numbers. The Ca is the space-averaged density, i.e. 

/ $.(r) = 0, (34) 

where Vu is the unit cell of the periodic structure, whose volume is denoted by Vu- In the 
case of periodic structures 

Cr^(ri + P, + P) = C™^(ri, r^) = ^^^(Arlr^) (35) 

where Ar = ri — r2 G i?^ and r2 G Vu- We introduce the inverse correlation function 
averaged over the unit cell by 

C„^(Ar) = ^ f (Ar|r2) (36) 
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with analogous definition for tlie correlation function Gap averaged over the unit cell (in 
terms of G™^). In Fourier representation from Eq. fl3T]) we have 25| 



C„,(k)G,^(k) = S^p- (37) 
We decompose -ff/^nc into two parts, 

Hfiucm, m = nam, m + ^nm, m. m 

The first term in the above equation is given by 

^g[C,0] = \ /^0a(k)C„^(k)0;3(-k) = ^ |^^,(k)a(k)7A.(-k), (39) 

where Cj(k) and are the eigenvalues and the eigenmodes respectively of the matrix C 

with the elements Capi}^)-, and summation convention for i is used. For brevity we introduced 
the notation = J j^^- In the next step we make an assumption that ATil^cM can be 
treated as a small perturbation. When such an assumption is valid, we obtain 
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28| 



PQ[Q ^ PQ,o[Q - log J D(l>^... J D0„e-^^« + {(3A'H)g + 0((/3AH)^). (40) 

where {■■■)g denotes averaging with the Gaussian Boltzmann factor e~^^^ . Eqs. (1301) - (HUj) 
allow for calculation of the fluctuation contribution to the grand potential in the lowest 
nontrivial order when the form of VL^o is known, and the form of Cap (Eqs. (l36p and (1271) ) is 
determined by a self-consistent solution of some approximate version of Eq. (!30|) . In general, 
each contribution to Eq. fHOj) depends on the mesoscopic length scale i?, but the i?-dependent 
contributions must cancel against each other to yield i?- independent Vt. By minimizing the 
density functional fHUl) we find the equilibrium structure. The main difficulty consists in 
determination of C*i(k). 

III. APPROXIMATE THEORY FOR WEAK ORDERING 

A. Grand potential and correlation functions in the case of weak ordering 

If ordering in the system occurs on a length scale larger than the size of particles, the 
local density approximation can be applied, and we assume 

mc}] = jj\ar\...Xn{v)\ (41) 
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where f'^{(i{r), ...,C„(r)) is the free-energy density of the hard-sphere system in which the 
volume fractions in the infinitesimal volume dr at r are Ci(r), Cn(i") • 
In this approximation we obtain the functionals 

/3^co[{C}] = 1 [ f ^V^"p{ri2)Uri)Ur2) - / /3A({C}) - / /3/x.C(r), (42) 



and 



/3/7^,„,[{C},0]= / /3[/.\(C(ri))-/ia.+ / ^.Ta.(^i2)C..(r2)]0«.(ri) (43) 

Jri '- Jr2 

+ 1 [ [ 0.(ri)Q(ri,r2)0^(r2) + E / ^^""-'"/^'^"^ V a.(r)-0.„(r). 

where V^"^^{ri2) is defined in Eq. flT^ and 

({C(r)}) = ^"^^'•»«'-)» (44) 

depends on the (local) composition of the mixture, but is independent of temperature. 
For inhomogeneous phases in soft matter systems (colloidal crystals for example), with 
position-dependent volume fractions, determination of the grand potential and the average 
distribution of particles within this theory is still very difficult. However, further simplifying 
assumptions can be made in the case of weak ordering. In the case of 'soft' crystalline phases 
with unit cells of the structure significantly larger than the size of the particles, particles 
(and the whole clusters) can fiuctuate around their average positions. The displacements 
of the particles from their average positions can be large, but the long-range order can be 
preserved. Averaging over such fiuctuations leads to smooth functions $Q(r) with small 
magnitudes $q, ^ (see Eq. (!32|) ). In the case of 'soft' crystalline phases the functional 
lOj) can be approximated by 

/3fi[{C° + Hr)}]/V ^ pn^iC' + Hr)}]/V (45) 



1=1 



In + Crik)G.ik) - 1 



8 

where Cf°{k) and Gi{k) denote the eigenvalues of the matrices C™(/c) and G{k) = C~^{k) 
respectively, with the element (a,/3) of the former given by 

C-,{k) = C%{k) + / (46) 

^ "'rev,, V,, 
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where + $a(r) characterizes local volume fraction of the a-th component in the inhomo- 
geneous ordered phase, and 



c'.,{k) = mm + /3/„^^({c"}). (47) 

In the disordered phase ^a(j) = for each component a, and C^^(A;) reduces to C°^(A;). 
Finally, 

G"^ = I G^p{k). (48) 

Recall that by construction of the mesoscopic theory on the length scale R, the cutoff ~ vr/i? 
is present in the integral in Eq. fHHj) . Recall also that from Eq.f l2ip and discussion in sec. 2.1 
it follows that J^Ga/3{k) = Gai3{0) differs from the microscopic correlation function at zero 
distance and is finite. Since the mesoscopic length scale is larger than the size of molecules 
and smaller than the length scale of the ordering but otherwise it is arbitrary, the above 
approximate version of the theory is valid as long as the i?-dependent terms in Eq. fHSl) are 
negligible compared to the dominant contribution. 

When ^ Cay Eq- (|42!) can be approximated by the expression 

/3fi,„[{C° + ^(r)}] = /3^]..[{C°}] + PnciiC' + <l>(r)}] (49) 

^„ ni Jr 

n>3 

with 

^^^g[{C' + Hr)}] = ^/^«l.(k)C'°,(A;)<l,(-k) = 1 ^',(k)C7°(fc)^,(-k), (50) 

where Ci{k) and ^j(k) are the eigenvalues and eigenvectors of the matrix C° respectively. 

In order to obtain approximation for Cap from Eq. fl30|) . we truncate the expansion of 
Hfiuc[{C}j {4>}] (sse Eqs. (!26l) and (143|) ) at the term 0{(f)^). Next, the four- and six-point 
correlation functions are approximated by products of two-point correlation functions, and 
after some algebra we obtain the approximate result, valid for periodic structures 

2C{k) = {&°{k) + A)[3I - G{k){&''{k) + A)] - P(A;) (51) 

where I is the unitary matrix (Jq,^ = S^^), and the (a, f3) element of the matrix A is 

= lPfap,M'}) j^G,,{k). (52) 
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Note that A is independent of k. Finally, the {a, (3) element of the matrix P is 



Rfh Pfh 

H / e Uo,^p^[r)Lxa^p^[r)Lxo,3P3[r). 



6 

When P is neglected in Eq.f lST]) . then we obtain a simple equation for Cik), analogous to 
the self-consistent Hartree approximation and Brazovskii theory 1^] generalized for mixtures 
(linear approximation with respect to /3/aj...Q,^), 

C{k) = &°{k) + A. (54) 

By inserting Eq. fl54l) into Eq. fl51l) one can easily check validity of the former when P is 
neglected. Note that the dependence on k in Eq. fl54|) is not changed compared to the 
MP result, because A is independent of k. However, when P is included in Eq.f l5ip . the 



dependence on k is different than in MP. In Ref.[25| the expansion of Hjiuc was truncated 
at the term 0(0^), which leads to less accurate approximation. However, at linear order 
in derivatives of the same self-consistent Hartree approximation (Eq. fl54p ) was obtained 
and used in further applications. 



B. Boundary of stability of the homogeneous phase in the Brazovskii-type ap- 
proximation 

The homogeneous system is unstable with respect to an infinitesimal fluctuation ipiif^) 
when Ci{k) < 0. The instability occurs when the fluctuation tpi{k) is excited and at the 
same time ipj{k) = for j ^ i. The boundary of stability of the homogeneous phase is given 
by 

detCik,) = l[Q{h) = 0, (55) 

i=l 

where kb corresponds to the highest temperature for which any instability occurs for given 
composition of the mixture {C°}- Since the temperature at the instability with respect to 
the concentration wave ipi{k) with the wave-number k, T{k), is given by det C{k) = 0, the 
(local) maximum condition dT/dk = is equivalent to 

ddetC(A;) 



dk 
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k=kh 



0. (56) 



If there are several solutions of Eqs. (155!) and (156|) . the one corresponding to the highest 
temperature for given composition {^^} determines the temperature and the wave number 
of the critical mode at the boundary of stability of the homogeneous phase with respect to 
concentration fluctuations with infinitesimal amplitudes. 

In MF Eg. (155 p reduces to detC°(/cb) = 0, which for fixed composition {(^} in an n- 
component system is an n-th order equation for (3 = 1/(A;bT) (see Eq. p7|) and note that 
/^/aio2 Ki^(^) independent of T). There are up to n solutions for T for fixed 
composition and k, depending on the interaction potentials. In the one-component system 
there is one such solution, and if it corresponds to k^ > 0, the corresponding line T{() is 



known as the A-line 



3, 
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30|. In multicomponent system we should speak about A-surface. 
In the mesoscopic theory the probability of a deviation from the average composition {(^} 
on the mesoscopic length scale, {0(r)}, is proportional to exp(— /3fico[{0(r) + C°}])- Inhomo- 
geneous distribution of the particles on the mesoscopic length scale can be more probable 
than the homogeneous states when (3Qco does not assume a minimum for {0(r) = 0}, i.e. 



when det C°(/c6) < 0. As discussed in Refs.|25l. l29L l3lj. at the A-line (or A-surface) a change 



from locally homogeneous to locally periodic structure occurs, because for det C^{kf)) < 
the waves with the wavelength 27i /kh (and infinitesimal amplitudes) are more probable than 
the constant volume fractions. However, in the presence of mesoscopic fluctuations the aver- 
aging over different waves of concentration may lead to position-independent average volume 
fractions, by which the stability of the disordered phase can be restored. 

The open question for a multicomponent system is whether solutions of Eqs. ( l55l) and ( 156|1 
exist beyond MF. To answer this question let us focus on C ( Eqs. ( 15^ . (H6l) and ( |52il ). and 
note that the matrix A is a linear combination of the integrals of the correlation functions, 

•*-<^' « 

det C{k) can be Taylor expanded near the minimum at k^, 

det C{k) ^ det C{kb) + C2{k - hf + 0{{k - k^f) (58) 



and from Refs. 17|, |25|, |32| we obtain 



n/R for kb = 

, + Oil R) for fcb > 

VdetC(fci,) 
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The cases kb = and > 0, corresponding to macro- and micro-phase separation respec- 
tively, are significantly different. In the case of separation into two homogeneous phases 
the temperature at the instability is shifted compared to the MF result, because the matrix 
elements of A are finite. Since A^i^ (x tc/R, the shift depends on the scale of coarse-graining. 
This approximation is oversimplified for precise determination of the spinodal line for the 
macroscopic phase separation. We can determine the upper bound of the shift, because 
R > o'ais- When the interaction potentials are such that Eg. (1561) leads to ki, > 0, then from 
Eg. (159!) it follows that the matrix A becomes singular if detC(/cb) — )■ 0, and Eg. (155!) has 
no solutions for finite C°(A;f,). Only for T — )■ 0, i.e. when C°(/cb) becomes singular as well, 
the solution of Eg. fl55p may exist. 

At first order transitions between inhomogeneous phases the grand-potential density takes 
the same value for different phases. Inhomogeneous phases may appear when the periodic 
mesoscopic densities are more probable than the uniform states, i.e. in the phase-space region 
where detC°(/cb) < 0. For such thermodynamic states the fluctuations on the mesoscopic 
length scale dominate, because the corresponding contribution to C is large enough to yield 
det C > despite det C*^ < 0. Thus, in order to obtain the phase diagram precisely, the 
fluctuation contribution to the grand potential cannot be neglected. We should note that 
the relatively simple MF stability analysis allows for the rough estimation of the phase space 
part where the structure may be inhomogeneous on the mesoscopic length scale. 



IV. TWO-COMPONENT MIXTURE 



The general results of the previous sections apply inparticular to two-component mix- 
tures. Two component mixtures were studied in Ref.jo] within the method of collective 
variables [Sj] under the assumption of hmomgeneous structure. Here we ara mainly inter- 
ested in inhomogeneous fluids. In the flrst step we need to determine the inverse correlation 
functions in MF approximation, and this is a subject of this section. The form of the free 
energy for hard spheres in two-component systems is known JbS], and was used recently 

351 ]. Hence, we can perform 



34, 



in the case of ionic systems with size asymmetry of ions 
more detailed analysis within the framework developed above, still for arbitrary form of the 
interaction potentials. 
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It is convenient to introduce 



= . (60) 



As a length unit we choose ai2, and the wave- numbers are in a^-z units. We shall use the 
index a = 1 for the larger, and the index a = 2 for the smaller particle. The asymmetry of 
the size of the particles can be characterized by 

n = ^ = 1 + 5 r2 = — = l-6, (61) 

0^12 Cri2 

where 0<6<1;6 = 0,1 for identical sizes of the particles and for point-like smaller 
particles respectively. We also introduce volume fraction of both components, 

C = Ci + C2. (62) 

We shall consider dimensionless correlation functions for local deviations of the volume 
fractions from the space-averaged values, $Q,(r) = Ca(i") ~ C«- 

A. Mean field approximation for arbitrary interaction potentials 

The partial inverse correlation functions in the disordered phase in MF approximation, 
C^^ik), are given in Eg. (1461) . Explicit expressions for dimensionless partial derivatives of 
the free-energy density for hard-sphere reference system, /a^(Ci)C2) = ^izP fapiCiXi) , are 
given in Appendix. In the case of ordering on the mesoscopic length scale we make the 
approximation ( |T71) for the interaction potentials in Eg. (1141) . When the microscopic structure 
is disregarded, the microscopic correlation function takes the form 

C^(ri-r2) = ^^(|ri-r2|-a,^), (63) 

and we approximate the interaction-potential density by 

K,iry)^olX1M . ('^]'Ml,^mEl_lA^, (54) 

Egs. (l55|) and (!56|) . from which the temperature at the boundary of stability of the disor- 
dered phase can be obtained, for the two-component mixtures take the forms 

P'^U{k) + I3k{k) + D = Q (65) 
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and 

pU'ik) + k'{k) = 0, (66) 
where the prime denotes a derivative with respect to k, and we have introduced 

D = det F U{k) = det Y*{k) (67) 

and 

K{k) = V*,{k)n, + V;,{k)n, - 2VUk)fu = DTr{V*{k)F-'). (68) 

By V*(A;) and F we denote the matrices with the (a,/3) element given by V*^{k) and f*p 
respectively. D can be directly calculated from F given in Appendix, and for any size ratio 
has the simple form 

CiC2(i-C)* 

Let us discuss conditions under which solutions of Eq. (j65l) for 1/(3, 

,,T,,.- f^'^-' , . (70) 



-K{k) ± ^JK{k)^ - ADU{k) 

are real positive numbers. Two types of interaction potentials can be distinguished: (i) 
U{k) < and (ii) U{k) > 0. 

(i) Interaction potentials such that U{k) < characterize, in particular, the Primitive 
Model of ionic systems (hard sphere and Coulomb potential) for k > [34]. Since D > 0, 
there is always one positive solution of Eg. (1651) . namely ^^72. 

(ii) The case of U{0) > corresponds, in particular, to stronger attraction (first moment 
of the interaction potential) between like particles than between particles of different kinds. 
For U{k) > the necessary condition for an instability with respect to a density wave with 
the wavelength k is K{k) < 0, because for K{k) > the solutions of Eq.( l65l) for (3, if exist, 
are negative. Another condition that must be satisfied for existence of solutions of Eq.( l65ll 
is K{kY — ADU{k) > 0. The higher temperature is for the solution A;bT2. 

In order to determine what kind of inhomogeneities appear in the system beyond the 
boundary of stability of the homogeneous phase we focus on the eigenvalues and eigenvectors 
of C^{k). The eigenvalues depend on the wave- number in a nontrivial way, 

^o(^) ^ CUk) + CUk)-stgn{CUk))B{k) ^^^^ 
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and 



with 



B{k) = JA^ + ACUky 



and 



Aik) = signiCUmC'22ik) - AV^)]- 
The corresponding eigenmodes have the forms 

^i(k) = S(A;)<li(k)-6(A;)<l2(k), 



^2(k) = 6(A;)<li(k) + S(A;)$2(k), 



where 



a{k) 



h{k) 



A{k) + B{k) 
2B{k) 

-A{k) + B {k) 
2B{k) 



1/2 



nl/2 



(72) 



(73) 



(74) 



(75) 
(76) 



(77) 



(78) 



The above expressions differ from the corresponding expressions derived for the primitive 
model of ionic systems in Ref. {34!, because here we consider volume fractions rather than 
number densities. Note that as the fluctuating field either the local number density or 
the local volume fraction can be chosen, and the physical properties like phase transitions 
should be independent of this choice. By changing Pa we should simultaneously rescale 
Capik) — > Caii{k)vaV(i, and the Eqs. (l55ll and (l56l) for the rescaled functions yield the same 
solutions. 

The eigenmodes represent two order-parameter (OP) fields, and in principle either one of 
them may lead to instability of the disordered phase for given (Ci, (,2)- For small-amplitude 
inhomogeneities the last term in Eq. f H9l) can be neglected, and we can limit ourselves to the 
Gaussian approximation in Eq. flSOp . The OP \E'i(kft) induces the instability when C?(/cfe) = 0. 
Recall that at the boundary of stability with respect to the fluctuation \E'i(kb) the other OP 
vanishes (it would yield a positive contribution to the grand potential). The requirement 
that ^^(kft) = for j 7^ i leads to the relation between the critical amplitudes (see Eqs. (176|) 
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and (175 



$i(kfe) 

Since -B(fc) ± ^{k) > 0, we have 



(79) 



si^n(<l2(k;,)<li(kfe)) 



-sign{CUh)) for ^^(A;;,) = 
[stgniCUh)) for ^"(A;,) = 



(80) 



Let us consider planar waves $j(r) = $.j(k;,) cos(k5 ■ r). The focal change of the volume 
fraction is C(r) - C° = ^lO^b) cos(kb ■ r) + ^2iK) cos(k;, ■ r). When ^2(^b)^i(^b) > 0, then 
dense regions of size 7c/kh where C(r) — C° > ^''^^ followed by dilute regions of similar size 
where C(r) - C < 0. When $2(kb)$i(k6) < 0, then regions of size rc/kh where the volume 
fraction of the first component is enhanced (Ci(r) > Ci) the same time the volume 

fraction of the second component is depleted (C2(r) < C2) followed by regions where 
Ci(r) < Cl and C2(r) > (2- Thus, when <l2(kfe)<li(kb) > global (for kb = 0) or local (for 
kh > 0) gas-liquid separation occurs, whereas the case ^2(}<-b)^i(Mb) < corresponds to 
global or local demixing. 



B. Examples 

For an illustration we consider here three simple examples of different types of interaction 
potentials, leading to different behavior. 

I. ^2*2 = = 0. 

In this case it follows from Eq. fl65|) that at the instability with respect to the k-mode the 
temperature is 

ksT = -YMIk^ (81) 

and the instability occurs only if Vii{k) < 0. The boundary of stability of the disordered 
phase corresponds to the minimum of V{i{k) (maximum of — (/c)). When V{i{r) < for 
all r, then the minimum of V^i{k) is assumed for A; = 0, because V^i* (0) = J^. V^^{r). However, 
for potentials that are positive for some distances and negative for another distances, like 
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. In the former 



the SALR potential, the minimum of Vii{k) may be assumed for k > 
case the macroscopic separation into phases rich- and poor in the first component occurs. 
In the latter case the instability is induced by the concentration wave of the first component 
with the wavelength 2ti /ki, - the regions of excess volume fraction are followed by regions of 
depleted volume fraction. The volume fraction of the second component is given in Eq.f l79|) . 
For the chosen interactions we have C'i2(^) > 0) ^^"^ critical mode is \^i(kfe) or \E^2(k6) 
for Cl^{k) + C'2°2(A;) > or C^^^k) + C^^{k) < respectively (see Eqs.([ZI])-([72D). Note that 
when Ciiikij) = and $2(^6) = 0, there is another instability, with respect to fluctuation 
of the first component only. The corresponding temperature, ksT = —Vii{kb)/fii, is lower 
than in Eq.(l8l]). 

II. 1/1*2 = and ^2*2 = 1/1*1. 

This case corresponds to the simplest interactions for which U{k) > 0. From Eq. fl65l) we 
obtain the instability with respect to the k-mode 

kT = -^^■■('^' (82) 
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/n + /l2- V(/n-/22)' + 4/: 

Again, the boundary of stability of the disordered phase corresponds to the minimum of 
1/1*1 (fc). Only for V^iik) < the instability with respect to the k mode can occur, because 
the denominator is positive. For attractive interactions the minimum of V{i{k) is for A; = 0. 
In this case Ci2{k) > too, and the critical mode is \I^i(kf,) or \I'2(kf,). In the first case 
<l'2(kb)$i(kh) < and in the second case $2(kb)$i(k;,) > 0, hence demixing and gas-liquid 
type separation occurs in the first and in the second case respectively. The sign of C*°(fc) 
depends on the volume fractions and T, hence gas-liquid separation and demixing occur for 
different parts of the phase diagram. 
III. 1^22 = V,\ = 0. 

In this case U{k) = —V{2{k) < and K{k) = — 2/1*2 1^*2 (A;). There is one positive and 
one negative solution of Eq. fl65l) for given k, and the positive solution takes the form 



for 1/1*2 (A;) < 



f'BT = { ^ff^'^^'' . . (83) 

-7^^ for y^ik) > 

It is instructive to consider a particular form of V*2{k). We choose square- well potential 



V*2{r) = -v9{r - l)9{a - r) 
25 



(84) 



-10 




FIG. 3: Interaction potential between particles of different kinds in Fourier representation, Eg. ([85 
for a = 2. The wave-number k is in crj^^ units, the potential is in v units. 



where a > 1 is the range of the potential and r is in ai2 units. In Fourier representation we 
have 

V* {k) 'i'^v{ak cos{ak) — k cos k + sink — sm{ak)) ^^^^ 

k^ 



V{2{k) is shown in Fig.3 for a = 2. For this potential the MF instabihty is given in Eq.( l83ll . 
in the upper or in the lower line for A; = or for k^ ~ 2.78 (the first maximum of V{2{k)) 
respectively. In the first case C°2(0) < and $2(0)$i(0) > 0, hence the gas-liquid separation 
occurs. In the second case Cl2{ki,) > and ^2{kb)^i{kh) < 0, hence local demixing occurs. 
In each case Cf2{kb)^2{kb)^i{—kh) < and the Gaussian part of the excess grand potential, 
Eq. (!50l) . can vanish (in this example Cfj(fc) > 0). The two surfaces in Eq.( !83l) are shown in 
Figs. 4 and 5 for equal sizes of particles, and in Fig.6 for 6 = 0.95 (o"i/(j2 = 39). The period 
of the most probable inhomogeneities is 2'K/kf) ^ 2.26cri2. For equal sizes ai2 = ci, but in 
the case of the large size asymmetry the characteristic extent of inhomogeneities is 1.16cri. 
This example is at the boundary of applicability of the mesoscopic description. 

For equal sizes gas-liquid separation occurs for low volume fractions, and for higher 
volume fractions the liquid phase undergoes periodic ordering. The MF instability with 
respect to periodic ordering signals tendency for formation of an "ionic" crystal, where 
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nearest-neighbors are of different kind 36|. For large size ratio the phase-separation is 
found for temperatures lower than the temperature at the A-surface for all volume fractions. 
Below the A-surface the system is inhomogeneous. Studies beyond MF are necessary to 
clarify whether transition between gas and lyotropic liquid crystal preempts the gas-liquid 
separation, or the transition between inhomogeneous fluids (containing clusters or other 
aggregates) occurs. 

Let us summarize the above examples. When the interaction potentials between like par- 
ticles are attractive for some distances and repulsive for different distances, and the Fourier 
transform assumes a negative minimum for A; > 0, then the system is inhomogeneous below 
the A-surface. Such behavior was found already in one-component systems. In addition, 
inhomogeneous structures can occur when particles of different kind attract each other and 
the Fourier transform of the interaction potential assumes positive maximum for A; > 0. 
Periodic ordering is enhanced when the size asymmetry increases, and the system becomes 
inhomogeneous even for low volume fractions for 5 > 0.9. Increasing tendency for clustering 
with increasing size asymmetry was observed in ionic systems, in mesoscopic theory 34] and 



in simulation studies 
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C. Beyond MF stability analysis 

In order to calculate the phase diagram and structure we need to calculate the correla- 
tion functions beyond MF. The explicit expressions for the inverse correlation functions in 
MF in principle allow for obtaining the correlation functions in the Brazovskii-type approx- 
imation by self-consistent solutions of Eq. (!5^ . In practice this is less trivial than in the 



one-component case 



32| 



The grand potential functional of the mesoscopic volume fractions can be obtained once 
the matrix C™(/i;) and its eigenvalues are determined from Eqs.f H5l) -(H9l). However, even at 
the MF level determination of the phase diagram is not quite trivial. The problem simplifies 
when we restrict our considerations to phases of particular symmetry. In one-component 



systems the MF phase diagram was determined for weak ordering in Ref . |25|, |39| , and in 
two-component case one can proceed in a similar way. The fluctuation contribution to the 
grand potential can be obtained from Eqs.f l45p -f lSU]) and f lM|) . We should stress that from 
sec. 3 it follows that the phase transitions between inhomogeneous phases are first-order. 
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FIG. 4: Surfaces representing dimensionless temperature ksT/v at the MF boundary of stability 
of the homogeneous phase for Vi2{k) given in Eq. (j84p with a = 2, and Vii{k) = V22ik) = 0, for 
equal sizes of particles. Instability with respect to gas- liquid separation is shown in top left panel 
and instability with respect to periodic ordering (A-surface) in top right panel. In bottom panel 
both surfaces are shown. Volume fractions of the two components (dimensionless) here are denoted 
by zl and z2 and in the text by Ci, C2 respectively. 

Temperature at the transition is shifted compared to the MF result, because the fluctuation 
contribution to the grand potential in different phases is different. 
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FIG. 5: Cross-sections through the surfaces shown in Fig. 3 (bottom panel). Left: for equal volume 
fractions of the two components, (i = C2- Right: for (2 = 0.35 — ("1. Solid and dashed lines represent 
gas-liquid separation and the A-surface (the upper and the lower line in Eq. (|83p ) respectively. 
kBT/v and the volume fraction are both dimensionless. 




FIG. 6: Surfaces representing dimensionless temperature ksT/v at the MF boundary of stability 
of the homogeneous phase for V*2{k) given in Eq. ([8^ with a = 2, and Vii{k) = V22(^) = 0, for 
S = 0.95 (o'i/cj2 = 39). The gas-liquid separation (lower surface) is metastable in MF for all values 
of the volume fractions. Volume fractions of the two components (dimensionless) here are denoted 
by zl and z2 and in the text by i^i, C2 respectively. 
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V. SUMMARY AND DISCUSSION 



We developed a mesoscopic DFT for inhomogeneous mixtures in terms of mesoscopic 
volume fractions. Local volume fraction is suitable for description of inhomogeneous dis- 
tribution of particles on mesoscopic length scales when particles of different species are of 
significantly different sizes. This is because derivation of a Landau-type theory from statis- 
tical mechanics requires expansions in local deviations of the volume fraction (or density) 
from the average value. Since volume fractions are less than unity regardless of the size 
of the particles, such expansions are justified from mathematical point of view. We per- 
formed precise coarse-graining by introducing the microscopic volume fraction, and then by 
averaging it over mesoscopic regions (see Eq.([H])) to obtain mesoscopic field. This field is 
next considered as a constraint on the microscopic states. Exact expression for the grand 
potential consisting of two parts is obtained. The first part contains contributions from the 
microscopic states that are compatible with the constraint imposed on the volume fractions 
on the mesoscopic scale - the mesoscopic volume fraction of each component must be equal 
to the ensemble average. This term has a form similar to the grand potential in standard 
DFT. The second term is the contribution resulting from mesoscopic fluctuations - that is, 
from microscopic states that are not compatible with the ensemble average of the volume 
fraction on the mesoscopic scale. This term in turn has a form similar to the Landau theory. 

In practice additional assumptions and approximations are necessary in order to obtain 
predictions for particular systems in the framework of this theory. Being interested in in- 
homogeneities on mesoscopic length scale, we assume that ordering occurs on length scales 
larger than the size of particles, and adopt local density approximation for the hard-sphere 
reference system. The functional can be further simplified under the assumption of weak 
ordering, by which we mean that local deviations of the volume fraction from the space 
averaged value are small. Under the two assumptions, often valid in soft matter, we obtain 
a functional (Eqs.( H5l) -( l50ll . ( 15^ and ( |52l) ) which is similar to an extension of the Landau 
theory for mixtures combined with simple DFT. The role of the OP's is played by linear 
combinations of the deviations of the local volume fractions from the space-averaged val- 
ues. Minima of the functional correspond to equilibrium structures. All parameters in the 
functional are expressed in terms of the free-energy density of hard-spheres and in terms of 
interaction potentials which can have arbitrary form. In the case of one-component systems 
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we obtain LGW or LB theory, when the interaction potential times the pair distribution 
function, V^", in Fourier representation is approximated by 



The LGW and LB theories correspond to fc;, = and kb > respectively. In original LB the- 
ory the second term in Eq.( l86l) is proportional to (fc^ — /c^)^, but for k ^ kb, i.e. for dominant 
wave-numbers, we have (/c^ — kf)'^ ^ {2kbY{k — kbY- From the theory developed in Ref . {25 1 
either the LGW or the LB theory is obtained as a further approximation, depending on the 
form of the interaction potential. We are not aware of extensions of the phenomenological 
LB theory to mixtures. From the theory developed in this work we can determine whether 
the system can separate in homogeneous phases, or whether inhomogeneous structures ap- 
pear on low-temperature side of the A-surface. This information can be obtained from the 
form of interaction potentials and from the size ratio of the particles, by performing stability 
analysis of the homogeneous phase (sees. 3 and 4). In Landau-type theories separation into 
homogeneous phases or periodic ordering is an apriori assumption. 

Despite strong assumptions and approximations, determination of the phase diagram in 
this theory is not easy, but we can draw some qualitative conclusions already from the 
relatively simple stability analysis. From the approximate form of the inverse correlation 
function (Eqs.( l5^ and ( 152|) ) it follows that instability with respect to periodic ordering 
[kb > 0) obtained in MF is removed by mesoscopic fluctuations, as was the case also in one 
component systems. The A-surface corresponding to the MF instability may be a border- 
line between homogeneous and inhomogeneous structure of the disordered phase. For high 
volume fractions formation of some kind of periodic crystal can be expected on this side 
of the A-surface which corresponds to inhomogeneous structure. For low volume fractions 
we may expect formation of ordered clusters, where domains with increased and depleted 
volume fraction of particular components are periodically ordered. However, the clusters 
('living polymers') can have different sizes and locations in space as a result of the meso- 
scopic fluctuations. Simple examples that illustrate the theory for two components indicate 
that formation of inhomogeneous distribution of particles is enhanced when the size ratio 
between the particles of different components increases. Further work is necessary for de- 
termination of phase diagrams for particular systems within the framework of the theory 
developed in this work, and for verification of validity of our approximations in different 
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experimental systems. It is important to note that formation (or not) of inhomogeneous 
distribution of particles in mixtures depends not only on the interaction potentials, but on 
a combined effect of interactions and size ratios. In future studies predictions of this the- 
ory for two components should be compared with predictions of the theory in which only 
the big particles are taken into account explicitly, and the small components lead only to 
solvent-mediated effective interactions. 
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VI. APPENDIX 



The partial derivatives of the dimensionless free-energy density for hard spheres, 
/a^(Ci)C2) = ^'"12/^/0/3 (Ci 5 C2) are obtained from the explicit expressions for the chemical 
potentials in the hard-sphere mixture derived in Ref.j^Sj. Direct differentiations lead to the 
following expressions 



and 



/aa(Cl5 C2) 



6 



TTVi 



l_ ^ _8_ ^ ISr^Xi + QrlX2 + r^X^ 



Ca 1-C 



;i-c)^ 



;i-c)= 



:i-c)^ 
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/r2(Ci,C2 



7r(rir2) 
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where 



and is defined in Eq.( l6T|) 
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